Magurano et at. Virology Journal 2012, 9:284 
http://www.virologyj.eom/content/9/1/284 



□ 



VIROLOGY JOURNAL 



RESEARCH 



Open Access 



Molecular epidemiology of measles virus in Italy, 
2002-2007 



Fabio Magurano , Claudia Fortuna, Antonella Marchi, Eleonora Benedetti, Paola Bucci, Melissa 
and Loredana Nicoletti 



len 



Abstract 

Background: The European Regional Office of the World Health Organization (WHO/Europe) developed a strategic 
approach to halt the indigenous transmission of measles in its 53 Member States by 201 5. In view of the goal of 
measles elimination, it is of great importance to assess the circulation of wild-type measles virus (MV). Genetic 
analysis is indispensable to understand the epidemiology of measles. 

Methods: Urine and saliva samples were collected between May 2002 and December 2007, in order to find the 
origins and routes of wild type measles virus circulation. RT-PCR was performed on a total of 414 clinical samples of 
patients from different Italian regions. The results confirmed the genome presence in 199 samples, out of which 
179 were sequenced. The sequences were genotyped by comparing the fragment coding for the carboxyl terminus 
of the nucleoprotein (450 nucleotides) with that one of the WHO reference strains. 

Results: From the year 2002 to the year 2007 phylogenetic analysis of measles sequences showed a predominant 
circulation of the D7 genotype in the Italian territory for the years 2002-2004. This genotype was replaced by D4 
and B3 genotypes in the biennium 2006-2007. During the same period C2, A, D5 and D8 genotypes were also 
detected. 

Conclusions: Genetic characterization of wild-type MV provides a means to study the transmission pathways of the 
virus, and is an essential component of laboratory-based surveillance. Knowledge of currently circulating measles 
virus genotype in Italy will help in monitoring the success of the measles elimination programme and will 
contribute to evaluate the effectiveness of future vaccination campaigns. 

Keywords: Measles outbreak, Molecular epidemiology, Genotype, Phylogenetic analysis 



Background 

Globally, measles morbidity and mortality rates have been 
dramatically reduced since 1963, as a result of the availabil- 
ity of a safe and effective vaccine and the implementation 
of enhanced vaccination strategies [1-3]. Interruption of in- 
digenous transmission of measles virus (MV) (Paramyxovir- 
idae; MorbOlivirus) has been reported for several countries 
[4]. Nonetheless, measles remains a leading cause of child- 
hood mortality worldwide, with an estimated 164,000 mea- 
sles deaths in 2008 (a 78% reduction compared to mortality 
rate in 2000), most of which took place in developing coun- 
tries, primarily because of underutilization of the vaccine 
[5,6]. Furthermore, large outbreaks continue to occur in 
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countries with high vaccination coverage, after importation 
of the virus from endemic regions [7] . WHO measles elim- 
ination plan in the European Region (EUR) is targeted by 
2015 [1]; thus, progress in measles control in this region is 
fundamental for the success of the measles elimination 
programme [8]. In 2002-2003, Italy experienced a large 
measles outbreak affecting mostly the Southern regions, 
with an estimated national incidence of 738/100,000 in 
2002 and 544/100,000 in 2003, in children below 15 years 
of age, corresponding to over 100,000 estimated cases in 
this age group [9-11]. The outbreak was due to poor 
measles vaccination coverage. In fact, in 2003 the Italian 
national rate of children vaccinated with one dose of 
measles-containing vaccine by 24 months of age was 77%, 
being significantly lower in Southern regions compared to 
Central and Northern areas [12,13]. 
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Measles is a statutory and notifiable disease in Italy. 
During a six-year period (2001-2006) an average of ap- 
proximately 5,400 cases was reported annually, with a 
range from 18,020 cases in 2002, to 215 cases in 2005. 
According to the National Elimination Plan, sensitivity, 
specificity, and timeliness of case reporting had to be 
improved. Therefore, in April 2007 an enhanced surveil- 
lance system was established [14], and a National Refer- 
ence Laboratory (NRL) was established at the Istituto 
Superiore di Sanita (ISS) in order to support cases ascer- 
tainment, confirm outbreaks/cases and determine the 
MV genotypes. 

Molecular epidemiology, i. e. genetic characterization of 
wild-type MVs combined with standard epidemiological 
methods, is an essential component of the laboratory- 
based surveillance. It is performed throughout the world 
by the WHO Measles and Rubella Laboratory Network, 
which serves 166 countries in all WHO regions. Virological 
surveillance has helped to document the interruption of 
transmission of endemic measles in some regions. It 
includes epidemiological investigation and laboratory con- 
firmation of all sporadic illnesses clinically consistent with 
measles. Moreover, laboratory-based surveillance for 



measles and rubella, including genetic characterization of 
wild-type viruses, permits to illustrate the progress towards 
measles elimination by differentiating viruses between indi- 
genous and imported. Molecular characterization of mea- 
sles viruses provides a valuable tool to measure the 
effectiveness of measles control programmes, and viro- 
logical surveillance needs to be expanded throughout the 
world, and conducted during all phases of measles control. 

Genetic characteristics of representative wild-type MVs, 
identified in Italy between 2003 and 2007, were analyzed 
in this study. 

Results 

Urine and saliva samples were collected between May 
2002 and December 2007 (Figure 1) in order to find ori- 
gins and routes of MV wild-type circulation. RT-PCR 
(reverse transcriptase polymerase chain reaction) was 
performed on a total of 414 clinical samples, 399 urine 
and 15 saliva respectively, collected from different Italian 
regions and coming from 414 patients. The results con- 
firmed the presence of the genome in 199 samples, out 
of which 179 were sequenced (Table 1). The sequences 
were genotyped by comparing the fragment coding for 
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Figure 1 Geographical distribution of measles virus (MV) genotypes isolated in Italy, between 2002 and 2004 (A), and between 2006 
and 2007 (B). 
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the carboxyl terminus of the nucleoprotein (450 nucleo- 
tides) with the one of the WHO reference strains. A rep- 
resentative set of sequences is listed in Table 2. 

The results of genetic analysis indicated that 21 out of 
23 specimens sequenced in 2002, either before or during 
the peak of the outbreak, belonged to genotype D7 and 
2 to genotype A (Figure 2). All the 21 sequences classi- 
fied as genotype D7 were closely related to each other, 
showing overall only a single nucleotide difference. This 
fact supports the hypothesis of a common origin of the 
epidemic. Furthermore, these strains showed a 100% iden- 
tity with those isolated in France in 2001 and 2003 (MVi/ 
Paris.FR/01/1, MVi/Lyon.FR/03/1), in Germany in 2000 
and 2001 (MVs/Mainz.DEU/06.00/1, MVi/Mainz.DEU/ 
07.01), in Canada in 2000 (MVs/Alberta.CAN/20.00/1), 
and in Belarus in 2003 (MVs/Minsk.BLR/17.03), and were 
99% identical to the MV D7 reference strain MVi/Illinois. 
USA/50.99. 

During 2003, 7 specimens were sequenced: all of them 
belonged to genotype D7. All these sequences were 
closely related to each other, showing 99-100% identity 
with the D7 strains circulating in Italy in 2002, and with 
the strains circulating in Europe during the years 2000- 
2003 (Figure 2). 

During 2004, 4 specimens were sequenced: 3 of them 
belonged to the genotype C2 and 1 to genotype A. MVs 
strains belonging to the genotypes C2 circulated in 
Germany in 1992, in Spain during 1992-1993, and in 
Great Britain during 1992-1995, suggesting a wide dis- 
tribution of this genotype throughout Europe. More 
recently, genotype C2 has been detected in Luxembourg 
as well as in Germany, the Czech Republic and Denmark; 
it was responsible for the epidemic that occurred in 
Morocco in 2003. The Italian C2 strain differed for a sin- 
gle nucleotide from that one isolated in Luxembourg dur- 
ing 1996-1997 [15], as well as from two strains imported 
into the USA in 1997-1998 [16], while it differed only for 
one nucleotide from the C2 Moroccan strain (Figure 2). 
No epidemiological links are available for these strains. 
The Italian C2 strains showed a sequence identical to the 



Moroccan and Spain strains circulating in 2003 (MVi/ 
Casablanca.MOR/19.03/2). 

All the 3 strains belonging to the genotype A, circulat- 
ing in 2002 and 2004, revealed a close relationship with 
the vaccine strains: the strains isolated in 2002 showed 
identical sequences and the other A strain showed the 
90% of identity with the vaccine sequence (Figure 2). 
The strains isolated in the period 2002-2004 came from 
Campania, Lazio, Umbria Marche and Veneto regions 
(Figure 1). 

No laboratory data are available for the year 2005 be- 
cause of the interruption of the surveillance activity that 
restarted in the year 2006. 

During 2006, 75 specimens were sequenced. The 
phylogenetic analysis showed that 59 sequences out of 
76 belonged to genotype D4 and 16 to the genotype B3 
(Figure 3). 

In 2007, 69 sequences were analyzed; 29 of them 
belonged to genotype D4, 35 to genotype B3, 4 to geno- 
type D5 and 1 to genotype D8. 

Therefore, the majority of MVs circulating between 
2006 and 2007 in Central and Northern Italy belonged 
to genotype D4 (n= 88), and appeared to be grouped 
into two different clusters: D4 1 (n= 61) and D4 2 (n= 27), 
respectively (Figure 3). 

The 73% (44 out of 60) of the D4 1 sequences were 
100% identical to each other; 22% (13 out of 60) differed 
from the latter for a single nucleotide; 3% (2 out of 60) 
showed overall a difference of three nucleotides, and 
1.6% of five nucleotides (1 out of 60). 

BLAST analysis in MeaNS database (Measles Nucleo- 
tide Surveillance; http://www.who-measles.org) showed 
that the D4 1 cluster sequences were for the 99% identi- 
cal to the strains that circulated in Europe and then 
imported into the USA in the years 2006 and 2007 [17]. 
The strains belonging to this cluster were first isolated 
in February 2006 and circulated until February 2007 in 
Trentino Alto Adige, Lazio [18] and Sardinia regions. 

The measles viruses belonging to the cluster D4 2 were 
first isolated in February 2006 and circulated throughout 



Table 1 Measles cases tested at the Italian National Reference Laboratory (ISS) during the period 2002-2007, and 
related genotypes 



N of cases Genotypes 



Year 


Tested samples 


RT-PCR positive 


Sequenced 


A 


B3 


C2 


D4 


D5 


D7 


D8 


2002 


38 


23 


23 


2 










21 




2003 


19 


/ 


/ 












7 




2004 


21 


4 


4 


1 




3 










2005 


no data 


no data 


no data 
















2006 


127 


93 


76 




16 




60 








2007 


209 


72 


69 




35 




29 


4 




1 


TOT 


411 


199 


179 


3 


51 


3 


89 


4 


28 


1 
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Table 2 List of the representative sequences analyzed between 2002 and 2007 



Year 


Name of strain 1 


N of sequences 


N of weeks 2 


Genotype 


GenBank accession number Total 


2002 


MVs/Roma.lTA/20.02/1 


1 


1 


D7 


JQ783051 23 




MVs/Roma.lTA/20.02/2 


3 


6 


D7 


JQ783052 




MVs/Ancona.lTA/26.02 


1 


1 


D7 


JQ783050 




MVs/Caserta.lTA/27.02/1 


5 


1 


D7 


JQ783049 




MVs/Caserta.lTA/27.02/2 


1 


! 


A 


JQ783038 




MVs/Ascoli.lTA/27.02/1 


8 




D7 


JQ783048 




MVs/Ascoli.lTA/30.02/5 


1 


! 


A 


JQ783039 




MVs/Ancona.lTA/45.02 


3 


1 


D7 


JQ783047 


2003 


MVs/CastelfrancoVeneto.lTA/09.03 


1 


! 


D7 


JQ783040 7 




MVs/Roma.lTA/09.03 


1 


I 


D7 


JQ783041 




MVs/Avezzano.lTA/14.03 


1 




D7 


JQ783042 




MVs/Roma.lTA/15.03 


1 




D7 


JQ783043 




MVs/Pistoia.lTA/23.03 


1 




D7 


JQ783044 




MVs/Pordenone.lTA/1 9.03/1 


1 


1 


D7 


JQ783045 




MVs/Pordenone.lTA/1 9.03/2 


1 


1 


D7 


JQ783046 


2004 


MVs/Todi.lTA/28.04/1 


3 


1 


C2 


EF490999.1 4 




MVs/Todi.lTA/28.04/4 


1 


1 


A 


EF469771.1 


2006 


MVs/Merano.lTA/33.06 


2 


, 


D4 1 


AM849094.1 75 




MVs/Roma.lTA/33.06/1 


38 


18 


D4 1 


EF533887.1 




MVs/Latina.lTA/35.06 


/ 


/ 


D4 1 


JQ783032 




MVs/Alghero.lTA/36.06 


3 


2 


D4 1 


JQ783033 




MVs/Roma.lTA/36.06 


6 


4 


D4 1 


JQ783034 




MVs/Roma.lTA/37.06/7 


1 


1 


D4 1 


JQ783035 




MVs/Pescara.lTA/39.06 


1 


! 


D4 1 


JQ783036 




MVs/Roma.lTA/43.06/2 


9 


10 


B3.1 


EF533886.1 




MVs/Roma.lTA/45.06/2 


1 


1 


D4 1 


JQ783037 




MVs/Pomezia.lTA/48.06/1 


2 




B3.1 


JQ783028 




MVs/Albano.lTA/48.06 


3 




B3.1 


JQ783029 




MVs/Aprilia.lTA/49.06 


1 




B3.1 


JQ783030 




MVs/Latina.lTA/50.06 


1 


! 


B3.1 


JQ783031 


2007 


MVs/Roma.lTA/01.07 


8 


14 


B3.1 


JQ783008 69 




MVs/Roma.lTA/02.07 


1 


1 


B3.1 


JQ783009 




MVs/Barletta.lTA/03.07 


1 


1 


B3.1 


JQ783010 




MVs/CapoD'Orlando.lTA/04.07 


1 


1 


B3.1 


JQ78301 1 




MVs/Latina.lTA/07.07 


1 


1 


B3.1 


JQ783012 




MVs/Salerno.lTA/07.07 


1 


! 


D4 1 


JQ782998 




MVs/Merano.lTA/07.07 


1 


1 


D4 1 


JQ782999 




MVs/Milano.lTA/09.07 


1 




D4 2 


JQ783007 




MVs/Roma.lTA/1 0.07/1 


5 




B3.1 


JQ783013 




MVs/Roma.lTA/1 0.07/3 


1 




B3.1 


JQ783014 




MVs/Grosseto.lTA/1 1 .07 


1 




B3.1 


JQ783015 




MVs/Parma.lTA/11.07 


/ 




B3.1 


JQ783016 




MVs/MilanoJTA/11.07 


1 




B3.1 


JQ783017 




MVs/Viterbo.lTA/1 1 .07 


1 




B3.1 


JQ783018 
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Table 2 List of the representative sequences analyzed between 2002 and 2007 (Continued) 



MVs/Roma.lTA/1 2.07 


1 


1 D8 


JQ783053 


MVs/Roma.lTA/1 3.07 


3 


1 D4 2 


JQ783006 


MVs/Ancona.lTA/1 3.07 


3 


4 B3.1 


JQ783019 


MVs/Viterbo.lTA/1 3.07 




1 B3.1 


JQ783020 


MVs/Ferrara.lTA/15.07 




1 B3.1 


JQ783021 


MVs/Bari.lTA/21.07 




1 B3.1 


JQ783022 


MVs/Padova.lTA/23.07 




1 B3.1 


JQ783023 


MVs/Bolzano.lTA/38.07 




1 D5 


JQ783027 


MVs/Lodi.lTA/38.07 




1 D5 


JQ783025 


MVs/Roma.lTA/38.07 




1 D5 


JQ783024 


MVs/Bologna.lTA/39.07 




1 D5 


JQ783026 


MVs/Cuneo.lTA/41.07 


12 


7 D4 2 


JQ783005 


MVs/MilanoJTA/41.07/2 




1 D4 2 


JQ783004 


MVs/Cuneo.lTA/46.07 




1 D4 2 


JQ783003 


MVs^rento.lTA/47.07/1 




1 D4 2 


JQ783002 


MVs/Trento.lTA/47.07/2 




1 D4 2 


JQ783001 


MVs/Torino.lTA/51.07/1 




1 D4 2 


JQ783000 



: The nomenctature indicates the geographical origin and the year in which each virus was first identified. Identical sequences collected in the same place were 
grouped, and a single representative sequence was submitted in GenBank. 



2 : Period of circulation of each strain from its first isolation. 

the entire year in Trentino Alto Adige, Lombardy, 
Piedmont, Lazio and Sardinia regions (Figure 1). These 
sequences showed a maximum difference of 2 nucleo- 
tides (<0.5%), difference that remained unchanged when 
the strains MVi/Kolar.IND/03.07/1 and MVs/London. 
GBR/25.07 were added to the alignment, while it showed 
an increase with the addition of the strain MVs/Gros- 
seto.ITA/a301 [19]. 

Instead, the 51 B3 strains, isolated during the biennium 
2006-2007, belonged to the subgroup B3.1 (Figure 3) [20]. 
The genetic diversity of all these B3.1 circulating strains 
was limited to a maximum difference of 4 nucleotides 
(<0.99%). In addition, the genetic analysis revealed that 
some of these Italian cases were also identical to a case 
from Albania in 2006 (MVi/Shkodra.ALB/44.06), and to a 
case from Great Britain in 2005 (MVs/London.GBR/ 
31.05). The sequences belonging to B3.1 genotype were 
first isolated in October 2006 and the relevant strains cir- 
culated until June 2007 in Veneto, Lombardy, Emilia 
Romagna, Tuscany, Marche, Lazio, Abruzzo, Apulia and 
Sicily regions. 

Discussion 

This study stands so far as the first description of the 
molecular epidemiology of MV in Italy. 

Before this study the information about strains circu- 
lating in Italy derived from some cases occurred in 
Italy and imported into the U.S.A., Great Britain and 
Luxembourg [21,22]. The genotypes described were D6 
(1996-97) and D8 (1999). 



In our study we found that 28 (82%) out of the 34 MV 
sequences detected from clinical specimens during 
2002-2004 belonged to genotype D7. Three (9%) were 
placed in genotype C2 and three (9%) were placed in 
genotype A (Figure 2). 

The results of our analysis indicate that genotype D7 
was the cause of the epidemic that occurred in 2002, 
and probably of the one occurred in 2003, with 18,020 
and 11,978 measles confirmed cases, respectively. There- 
fore, genotype D7 was the endemic MV in the years 2002 
and 2003. The sequence analysis of the Italian cases 
showed that all the strains identified were strongly related 
to each other, and closely related to the European strains. 

Genotype D7 circulated in Great Britain and in Australia 
during the 1980s [23,24]. Chains of transmission of this 
genotype have been associated, through epidemiological 
investigations, to Sweden and other European countries 
[25]. In the early 2000s, this genotype replaced the geno- 
types C2 and D6, becoming the most commonly isolated 
genotype in Germany and France [26,27]. This shift was 
demonstrated in Germany during 2000-2001 [28]. Geno- 
types shift takes place in those countries where the inter- 
ruption of endemic transmission occurs only for short 
periods, due to measles sub-optimal control programmes. 
Indeed, failure in maintaining high levels of population 
immunity results in the accumulation of susceptible indivi- 
duals, thus creating favourable conditions for a rapid trans- 
mission of a newly introduced genotype. 

The comparison of the Italian A genotype MVs with 
those deposited in Means database revealed a close 
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MVs/California. USA/29.02 EU1 39095. 1 

- MVs/CastelfrancoVeneto. ITA/09.03 
MVi/Mainz.DEU/07.01 AF474938.1 
MVs/Roma. ITA/09.03 
MVs/Caserta.lTA/27.02/1 (n. 5) 
MV/Marseille.FR/03/4 AY585734.1 
MVs/Ancona.lTA/26.02 
MVs/Ancona.lTA/45.02 (n. 3) 
MVs/Roma.lTA/15.03 

- MVs/Pordenone.lTA/1 9.03/2 
MVi/Paris.FR/01/1 AY585732.1 
MV/Marseille.FR/03/2 AY585731.1 

- MVs/Roma.lTA/20.02/1 
MVs/Pordenone.lTA/19.03/1 
MVs/Mainz.DEU/06.00/1 AF277805.2 
MVs/Avezzano.lTA/14.03 
MVs/Pistoia.lTA/23.03 
MVs/Roma. ITA/20.02/2 (n.3) 
MVs/Ascoli.lTA/27.02/1 (n. 8) 

MVs/Alberta.CAN/20.00/1 AF41 1001.1 
MVi/Lyon . FR/03/1 A Y585729. 1 

- MVs/Minsk.BLR/17.03 AY523580.1 
MVi/lllinois. USA/50.99 AY037020.1 [D7] 

Schwarz-vaccine U03668. 1 

MVs/Caserta.lTA/27.02/2 

MVs/Ascoli.lTA/30.02/5 

MVs/Todi.lTA/28.04/4 

Edmonston-wt.USA/54 U01987.1 [A] 
Zagreb-vaccine AY486084.1 

- Maryland. USA/77 M89921 [C2] 
MVs/Barcelona.SPA/42.03 EU627192. 1 
MVs/Todi.lTA/28.04/1 (n.3) 
MVs.Taroudant.MOR/1 9.04/3 DQ779212.1 
MVi/Sefrou.MOR/06.03 DQ267525.1 
MVi/Casablanca.MOR/1 9.03/2 DQ267528.1 

- MVi/Reuler.LUX/16.96/1 AJ244032.1 
MVi/lndiana.USA/16.98 AY037011.1 



D7 



A 



C2 



Figure 2 Phylogenetic analysis according to the partial sequence of the N gene of measles virus (MV) strains identified in Italy, from 
2002 to 2004. Numbers in parenthesis indicate the number of additional identical strains from the same region. Significant bootstrap values 
(>80%) are indicated. 
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MWFIonda.USA/36.05 EU139106.1 
-MVs/Glasgow.GBR/20.07 EF653361.1 
MVs/Dublin.iRL/13.06 EF607929 1 
MVi/Lyon.FRA/24.06/2 EF428178.1 
MVs/Novi Sad.SRB/05.07 AM849060.1 
MViM/ashingtori.USA/46.05 DQ888752. 1 
MWBucharest.ROU/47.05/1 AM849092. 1 
MVs/Latina.lTA/35.06 (ri = 7) 
MVs/Roma.lTA/36.06 <n = 6) 
MVs/Sarajevo.BIH/1 1 .07/1 AM849061. 1 
— MVs/Merano.lTA/33.0G (n = 2) 
|MVs/Ber1in.DEU/16.05 AM849083.1 
-lMVs/Aarau.SWl/15.05 AM849090.1 
'MVs/Guimaraes.POR/29.05/1 AM849058.1 
MVs/Rom a. ITA/33.06/1 (n = 38) 
MVs/Salerno.1TA/07.D7 
MVs/Alghero.lTA/36.0G (n = 3) 
MVs/Pescara.lTA/39.06 
MVs/Merano.lTA/07.07 
MVs/Tarragona.ESP/1 0.07/1 FJ938034. 1 
MVi/Barcelona. ESP/2 07/2 FJ938033. 1 
MVs/Roma.lTA/37.06/7 
MWWashington USA/49.06 EU139111. 1 

MVs/Roma.lTA/45.06/2 

— MVs/Montreal.CAN/89 U01976.1 [D4] 

mVs/Grosseto ITA/a301 EF017348.1 



|— MVs/Kangra.lnd/3906/2 FJ223148.1 
MVs/Cuneo.lTA/46.07 (n = 7) 
MVs/Milano.lTA/41.07/2 
MVs/London.GBR/25 07 EU585740.1 
MVsfTorino.lTA/51.07/1 
MVs/Cuneo.lTA/41.07 fn = 12) 
MVs/Emmendingen.DEU/31.07 HQ704316 
MVs/Milano.lTA/09.07 
MVs/Harrow GBR/26.07 EU585741.1 
MVs/Roma.lTA/13.07 (n = 3) 
MVsfTrento.lTA/47.07/1 

MVs/Trento.lTA/47.07/2 

MVs/Kolar.lND/44.06 EU812271.1 
|-MWHetauda.NEP/2 99 AJ250071 1 [D8] 

MVs/Mosccw.RUS/33.07 EU597264 
MVs/Bagalkot.lND/48.06 EU812272.1 
MVs/Roma.lTA/12.07 
MVs/London GBR/10.07 EF554312.1 
■MVi/Bangkok THA/12 93 AF079555.1 [D5] 
MVs/Roma.lTA/38.07 
MVs/Lucerne.CHE/16.07 M'15615 

MVs/Lodi/ITA/38.07 

98 MVs/Bologna.lTA/39.07 

MVs/Osaka.JPN/31.07/1 AB593391 1 
MVs/Pensylvania USA/34.07 EU37S754 
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Figure 3 Phylogenetic analysis according to the partial 
sequence of the N gene of measles virus (MV) strains identified 
in Italy, in the biennium 2006-2007. Numbers in parenthesis 
indicate the number of additional identical strains from the same 
region. Significant bootstrap values (>80%) are indicated. M*: MeaNS 
id number. 



relationship with the Schwarz Italian vaccine strains, as 
well as with the wild-type genotype A strain (Figure 2). 
Unfortunately, no information about the vaccination sta- 
tus of these patients was available. However, the close 
genetic relationship with the vaccine strains indicates 
that these cases should be caused by vaccine virus. 

Phylogenetic analysis for the sequences of the Italian 
cases in the years 2006-2007 revealed a co-circulation of 
D4 and B3 genotypes. In addition, we identified an out- 
break caused by genotype D5, and only a single case 
belonging to D8 genotype in 2007 (Table 2; Figure 3). 

Eighty-nine (61%) out of the 145 MV strains detected in 
clinical specimens, in the years 2006-2007, belonged to 
genotype D4. This genotype was endemic on the Indian 
subcontinent, as well as in East and South Africa [25,29]. 
Genotype D4 has been repeatedly identified in the WHO 
Eastern Mediterranean Region [30], as well as in out- 
breaks and sporadic cases in several European countries, 
including Germany, Turkey, Spain, Great Britain, Croatia 
and Russia [17,31,32]. 

Italian D4 strains grouped into two different clusters 
of common origin (Figure 3). The strains belonging to 
cluster D4 1 were isolated from August 2006 to February 
2007, and were closely related to the strains circulating 
in Europe and USA from late 2005 to early 2007. 
According to the analysis made by Kremer et al. [17] 
these strains belonged to a group described as 'European 
group 1'. Moreover, this group was responsible for a 
large outbreak in Romania, which included >8,000 cases 
and lasted from December 2004 until early 2007 [17]. 

D4 Italian strains belonging to cluster D4 2 were isolated 
from November 2006 to December 2007, and were closely 
related to the Grosseto strain [19] and to the strains circu- 
lating in India, Great Britain and Germany in the same 
period. D4 2 was responsible for the outbreak occurred in 
Piedmont in 2007 having been imported from the United 
Kingdom [14] and differed from the European groups that 
circulated in the biennium 2005-2006 [17]. 

In 2007 the strain D5 caused a small outbreak in Italy, 
and the phylogenetic analysis suggests that it was prob- 
ably imported from Japan. In fact, despite unavailability 
of epidemiological links, this strain caused an outbreak 
in 2007 in Hokkaido district (Japan) [33]. 

In 2007 we found a single case belonging to the D8 
genotype (Figure 3), whose strain has been isolated in 
Russian [34] and was endemic in India in the same period 
[35] but the source of the virus remained unknown. 

From October 2006 to December 2007 the circulation 
of B3 strain, the same strain that caused the Apulia out- 
break in 2006 [36], was observed with the first case iden- 
tified. The strain was identified for the first time in the 
Lazio Italian region. This could suggest the place where 
the B3 outbreak started, but no epidemiological data 
supported this hypothesis. 
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Fifty-one (35%) of the MV strains detected in clinical 
specimens in the years 2006-2007 belonged to genotype 
B3. All these strains belonged to genotype B3.1 and were 
closely related to the strains isolated in Albania, Switzerland 
and circulating in Great Britain in the same period, suggest- 
ing a possible route of introduction by the latter (Figure 3) 
[14] although the origin of B3 outbreak was not known. 

In summary, according to our analyses we can assert 
that D7 strain was endemic in Italy from 2002 to 2004. 
In 2006 it was replaced by D4 and B3 strains that circu- 
lated also in 2007. 

In conclusion, a continuous improvement of labora- 
tory surveillance is needed, particularly of the molecular 
laboratory. This is significantly important in countries 
where the measles elimination phase, such as Italy, is tak- 
ing place and increasing efforts to obtain appropriate spe- 
cimens from each chain of transmission are necessary. 

Materials and methods 

Measles case definitions 

Measles cases were clinically defined i.e., generalized rash, 
lasting at least 3 days, accompanied by fever > 38°C, and 
cough, coryza or conjunctivitis. 

Laboratory confirmation was attained by the determin- 
ation of measles specific IgM antibodies through enzyme 
immunoassay, and/or by the detection of measles virus 
genome in saliva or urine by RT and hemi-nested PCR. 

Specimens 

Oral-fluid specimens were obtained by salivary swabs 
(Salivette, Sarstedt Company). Saliva was collected after 
centrifugation at 500 g for 10 minutes, and then stored 
until tested. 

Urine samples were collected within 7 days from rash 
onset. Urine sediment was obtained after centrifugation at 
500 g for 10 minutes, washed two times in sterile PBS, and 
resuspended in a final volume of 0.5-1 ml. Sediments were 
stored at -80°C. The tests on urine and salivary samples 
were performed under the indications of the "National 
Plan of Elimination of Measles and Congenital Rubella". 
The Plan has been elaborated by the Ministry of Health 
according with the WHO indications, and approved by all 
Regional Health Authorities. Our laboratory belongs to the 
WHO European Regional Network of National Measles 
and Rubella Reference Laboratories. Goal for this network 
is the measles elimination in Europe in 2015. Before col- 
lecting samples, subjects, or their parents in case of chil- 
dren, had to sign the "informed consensus". No formal 
approval from an ethics committee has been required; 
however the study complies with the Helsinki declaration. 

RT-PCR amplification and sequencing 

Total RNA was extracted using QIAmp Viral RNA Mini 
Kit (Qiagen) for saliva samples and RNeasy mini Kit 



(Qiagen) for urine samples, as per manufacturer proto- 
cols. Nucleic acid was tested by RT-PCR using a hemi- 
nested protocol [24] directed to a highly conserved part 
of the MV RNA, which is located on the N gene. Kit 
Superscript One-Step RT-PCR kit with Platinum Taq 
(Invitrogen) was used for RT-PCR reaction. 

Nucleotide Sequence Analysis 

Both strands of amplified products were sequenced by 
Macrogen Inc. (Seoul, Korea), using MVF2 and MVB1 
primers. 

Genomic sequences of reference strains used for gen- 
etic analysis of wild-type measles viruses were obtained 
from database by accession number [37,38]. Nucleotide 
sequences were aligned with the CLUSTAL W (BioEdit) 
software [39]. Phylogenetic trees were constructed using the 
nucleotide Kimura-2 parameter and the neighbour-joining 
method. Bootstrap analyses were performed through 1,000 
resampling of the data sets. The neighbour-joining method 
[40] was implemented by using MEGA-5 [41]. 

The sequences had been submitted to GenBank, and 
the corresponding accession numbers are given in Table 1. 
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